clear global;clc;
x1=-3;
x2=2;
y1=5;
po1=[x1,0];
po2=[x2,0];
po3=[0,y1];

tmp=[po1;po2;po3;po1].';
plot(tmp(1,:),tmp(2,:))

% xlist=linspace(x1,x2);
% for ni=1:length(xlist)
% ylist(ni) = ymax1(xlist(ni),x1,x2,y1);
% end
% 
% figure(2)
% plot(xlist,ylist);

testfun=@(x,y) 8*x.*y+x.^2+y.^2;

integral2(testfun,x1,0,0,@(x) x.*abs(y1./x1)+y1)+...
integral2(testfun,0,x2,0,@(x) x.*(-abs(y1./x2))+y1)

% delta_e=1;
% q2=delta_e*ymy_gauss_int(testfun,4)








